|
|
Computational Modelling of Self- oscillations of the Human Vocal Folds
Hájek, Petr ; Šidlof,, Petr (referee) ; Radolf, Vojtěch (referee) ; Švancara, Pavel (advisor)
The presented dissertation thesis deals with a simulation of the human phonation in terms of latest theories. Phonation is considered here as a bi-directional fluid-structure-acoustic interaction, where the interaction between all three physical domains occurs due to the unsteady viscous compressible Navier-Stokes equations. There is a solid knowledge background in the first part of the thesis. It concerns the latest concepts in computational modeling of the human phonation, the most important and recent theories about the human voice production and some key aspects of the human anatomy, physiology and pathology. Also voice assessment is discussed. The second part of the thesis describes an in-depth analysis of a phonation simulation in a planar computational model. The basic concepts proceed from algorithms developed in the Institute of Solid Mechanics, Mechatronics and Biomechanics. Created models are able to reproduce sounds of all Czech vowels and the most common evaluated parameters very close to physiological ranges. The simulated pathology, Reinke's edema, is demonstrated in order to explore its influence on the vowel sound. The third part focuses on modeling of phonation in a spatial computational model. All Czech vowels are simulated also here and compared to the planar model and to actual measurement. The spatial model serves as the starting point to modeling of a longitudinal pretension incorporated in the vocal folds. In the last part of the thesis, a modeling of the phonation with vocal folds pretension is investigated. Although the models are tuned to a rather soft phonation, the results are in agreement with the relevant physiologic phenomena. While the spatial model is highly computationally expensive, a hybrid planar model with pretension is proposed. A special attention is paid to the analysis of self-sustained oscillation of the vocal folds. It is shown, the planar model cannot reproduce such kind of oscillation in the actual version, albeit time of oscillation was considerably extended. On the other hand, oscillation of the spatial vocal folds are stabilized without effects accompanying subduing of oscillation. It can be supposed that the spatial model is able to reproduce self-sustained oscillation as a basic principle present during the human phonation.
|
|
|
Compensatory Vocal Folds for Source Voice Generation: Computational Modeling of Vocal Folds Function
Matug, Michal ; Vampola, Tomáš (referee) ; Horáček, Jaromír (referee) ; Švancara, Pavel (advisor)
This doctoral thesis focuses on computational modelling of human vocal folds and vocal tract functions using finite element method (FEM). Human voice is crucial in human communication. Therefore one of the main targets of current medicine is creation of artificial vocal folds, which would substitute the original vocal folds. The computational modelling can be used to understand principles of voice production, determination of parameters that the artificial vocal folds have to meet and verification of their functionality. First part of this thesis focuses on modelling of human voice creation by whisper. Influence of intraglottal gap on eigenvalues distribution for individual vowels was analysed using FEM vocal tract and trachea model. Further there is presented two-dimensional (2D) finite element model of the flow-induced self-oscillation of the human vocal folds in interaction with acoustic spaces of the vocal tract. The 2D vocal tract model was created on the basis of converting the data from magnetic resonance images (MRI). Explicit coupling scheme with separated solvers for structure and fluid domain was used for modelling of the fluid-structure interaction. Created computational model comprises: large deformations of the vocal folds tissue, contact between vocal folds, fluid-structure interaction, morphing the fluid mesh according to the vocal-fold motion (Arbitrary Lagrangian-Eulerian approach), unsteady viscous compressible or incompressible airflow described by the Navier-Stokes equations and airflow separation during glottis closure. This model is used to analyse the influence of stiffness and damping changes in individual vocal fold tissue layers (in particular in superficial lamina propria). Part of this computational analysis is also comparison of vocal folds behaviour for compressible and incompressible flow model. Videokymograms (VKG) are subsequently created from obtained results of FEM calculations which enable to compare individual variants between themselves and with motion of real human vocal folds. In next part of this thesis is presented three-dimensional (3D) finite element model of the flow-induced self-oscillation of the human vocal folds. This 3D model was created from a previous 2D model by extrude to the third direction. Using this model was again compared influence of compressible and incompressible flow model on vocal folds motion and generated sound by using videokymograms and acoustic spectra. The last part of this thesis focuses on the possibility to replace missing natural source voice in form reed-based element. Behaviour of reed-based element was analysed using computational modelling and using measurements on experimental physical model. The physical model enables changes in setting gap between reed and reed stop and performing acoustical and optical measurements.
|
|
|
Computational Modelling of Self- oscillations of the Human Vocal Folds
Hájek, Petr ; Šidlof,, Petr (referee) ; Radolf, Vojtěch (referee) ; Švancara, Pavel (advisor)
The presented dissertation thesis deals with a simulation of the human phonation in terms of latest theories. Phonation is considered here as a bi-directional fluid-structure-acoustic interaction, where the interaction between all three physical domains occurs due to the unsteady viscous compressible Navier-Stokes equations. There is a solid knowledge background in the first part of the thesis. It concerns the latest concepts in computational modeling of the human phonation, the most important and recent theories about the human voice production and some key aspects of the human anatomy, physiology and pathology. Also voice assessment is discussed. The second part of the thesis describes an in-depth analysis of a phonation simulation in a planar computational model. The basic concepts proceed from algorithms developed in the Institute of Solid Mechanics, Mechatronics and Biomechanics. Created models are able to reproduce sounds of all Czech vowels and the most common evaluated parameters very close to physiological ranges. The simulated pathology, Reinke's edema, is demonstrated in order to explore its influence on the vowel sound. The third part focuses on modeling of phonation in a spatial computational model. All Czech vowels are simulated also here and compared to the planar model and to actual measurement. The spatial model serves as the starting point to modeling of a longitudinal pretension incorporated in the vocal folds. In the last part of the thesis, a modeling of the phonation with vocal folds pretension is investigated. Although the models are tuned to a rather soft phonation, the results are in agreement with the relevant physiologic phenomena. While the spatial model is highly computationally expensive, a hybrid planar model with pretension is proposed. A special attention is paid to the analysis of self-sustained oscillation of the vocal folds. It is shown, the planar model cannot reproduce such kind of oscillation in the actual version, albeit time of oscillation was considerably extended. On the other hand, oscillation of the spatial vocal folds are stabilized without effects accompanying subduing of oscillation. It can be supposed that the spatial model is able to reproduce self-sustained oscillation as a basic principle present during the human phonation.
|
|
|
Three-dimensional numerical analysis of Czech vowel production
Hájek, P. ; Švancara, P. ; Horáček, Jaromír ; Švec, J. G.
Spatial air pressures generated in human vocal tract by vibrating vocal folds present sound sources of vowel production. This paper simulates phonation phenomena by using fluid-structure-acoustic scheme in a three-dimensional (3D) finite element model of a Czech vowel [o:]. The computational model was composed of four-layered M5-shaped vocal folds together with an idealized trachea and vocal tract. Spatial fluid flow in the trachea and in the vocal tract was obtained by unsteady viscous compressible Navier-Stokes equations. The oscillating vocal folds were modelled by a momentum equation. Large deformations were allowed. Transient analysis was performed based on separate structure and fluid solvers, which were exchanging loads acting on the vocal folds boundaries in each time iteration. The deformation of the fluid mesh during the vocal fold oscillation was realized by the arbitrary Lagrangian-Eulerian approach and by interpolation of fluid results on the deformed fluid mesh. Preliminary results show vibration characteristics of the vocal folds, which correspond to those obtained from human phonation at higher pitch. The vocal folds were self-oscillating at a reasonable frequency of 180 Hz. The vocal tract eigenfrequencies were in the ranges of the formant frequencies of Czech vowel [o:] measured on humans, during self-oscillations the formants shifted to lower frequencies.
|
|
|
Influence of tissue changes in superficial lamina propria on production of Czech vowels
Hájek, P. ; Švancara, P. ; Horáček, Jaromír ; Švec, J.
Superficial lamina propria (SLP) is a water-like vocal fold (VF) layer located directly under overlying epithelium. Its material properties affect VF motion and thus resulting spectrum of produced sound. Influence of stiffness and damping of the SLP on sound spectrum of Czech vowels is examined using a two-dimensional (2D) finite element (FE) model of a human phonation system. The model consists of the VF (structure model) connected with an idealized trachea and vocal tract (VT) (fluid models). Five VTs for all Czech vowels [a:], [e:], [i:], [o:] and [u:] were used and their geometry were based on MRI data. Fluid flow in the trachea and VT was modelled by unsteady viscous compressible Navier-Stokes equations. Such a formulation enabled numerical simulation of a fluid-structure-acoustic interaction (FSAI). Self-sustained oscillations of the VF were described by a momentum equation including large deformations and a homogeneous linear elastic model of material was used. Fluid and structure solvers exchange displacements and boundary forces in each iteration. During closed phase VFs are in contact and fluid flow is separated. We can observe that both the damping and the stiffness of the SLP substantially influence the amplitude and frequency of VFs vibration as well as the open time of the glottis.\n
|
|
|
Effect of turbulence in FE model of human vocal folds self-oscillation
Hájek, P. ; Švancara, Pavel ; Horáček, Jaromír ; Švec, J.G.
The purpose of the study is to determine whether a turbulence model in fluid flow calculation affects the vocal folds (VF) vibration and the acoustics of human vocal tract (VT). The objective is examined using a two-dimensional (2D) finite element (FE) model of the fluid-structure-acoustic interaction for self-sustained oscillations of the VF. The FE model consists of the models of the VF, the trachea and a simplified model of the human VT. The developed FE model includes large deformations of the VF tissue and VF contact interrupting the airflow during glottis closure. The airflow is modelled by the unsteady viscous compressible Navier-Stokes equations, without and with the Shear Stress Transport (SST) turbulence model. Fluid-structure interaction (FSI) and morphing of the fluid mesh are realized using Arbitrary Lagrangian-Eulerian (ALE) approach. The method is applied on the FE model of the VT shaped for the Czech vowel [a:]. Also effect of varying stiffness of the superficial lamina propria (SLP) is analyzed. The numerical simulations showed that considering of the turbulence affects mainly higher frequencies apparent in a frequency spectrum of the VT acoustics.
|
|
|
Numerical simulation of the effect of stiffness of lamina propria on the self-sustained oscillation of the vocal folds
Hájek, P. ; Švancara, P. ; Horáček, Jaromír ; Švec, J. G.
A two-dimensional (2D) finite element (FE) model of the fluid-structure-acoustic interaction during selfsustained oscillation of the human vocal folds (VF) is presented in this paper. The aim is to analyze the effect of stiffness of lamina propria on VF vibrations. Such stiffness change can be caused by some VF pathologies. The developed FE model consists of the FE models of the VF, trachea and a simplified human vocal tract. The vocal tract model shaped for simulation of phonation of Czech vowel [a:] was created by converting data from the magnetic resonance images (MRI). The developed FE model includes VF contact, large deformations of the VF tissue, fluid-structure interaction (FSI), moving boundary of the fluid mesh (Arbitrary Lagrangian-Eulerian (ALE) approach), airflow separation during the glottis closure and solution of unsteady viscous compressible airflow described by the Navier-Stokes equations. The numerical simulations showed that higher values of lamina propria Young's modulus (stiffer lamina propria) result in a decrease of the maximum glottis opening. Stiffer lamina propria also requires the use of higher subglottal pressure to initiate self-sustained vibration of the VF.
|
|
|
Compensatory Vocal Folds for Source Voice Generation: Computational Modeling of Vocal Folds Function
Matug, Michal ; Vampola, Tomáš (referee) ; Horáček, Jaromír (referee) ; Švancara, Pavel (advisor)
This doctoral thesis focuses on computational modelling of human vocal folds and vocal tract functions using finite element method (FEM). Human voice is crucial in human communication. Therefore one of the main targets of current medicine is creation of artificial vocal folds, which would substitute the original vocal folds. The computational modelling can be used to understand principles of voice production, determination of parameters that the artificial vocal folds have to meet and verification of their functionality. First part of this thesis focuses on modelling of human voice creation by whisper. Influence of intraglottal gap on eigenvalues distribution for individual vowels was analysed using FEM vocal tract and trachea model. Further there is presented two-dimensional (2D) finite element model of the flow-induced self-oscillation of the human vocal folds in interaction with acoustic spaces of the vocal tract. The 2D vocal tract model was created on the basis of converting the data from magnetic resonance images (MRI). Explicit coupling scheme with separated solvers for structure and fluid domain was used for modelling of the fluid-structure interaction. Created computational model comprises: large deformations of the vocal folds tissue, contact between vocal folds, fluid-structure interaction, morphing the fluid mesh according to the vocal-fold motion (Arbitrary Lagrangian-Eulerian approach), unsteady viscous compressible or incompressible airflow described by the Navier-Stokes equations and airflow separation during glottis closure. This model is used to analyse the influence of stiffness and damping changes in individual vocal fold tissue layers (in particular in superficial lamina propria). Part of this computational analysis is also comparison of vocal folds behaviour for compressible and incompressible flow model. Videokymograms (VKG) are subsequently created from obtained results of FEM calculations which enable to compare individual variants between themselves and with motion of real human vocal folds. In next part of this thesis is presented three-dimensional (3D) finite element model of the flow-induced self-oscillation of the human vocal folds. This 3D model was created from a previous 2D model by extrude to the third direction. Using this model was again compared influence of compressible and incompressible flow model on vocal folds motion and generated sound by using videokymograms and acoustic spectra. The last part of this thesis focuses on the possibility to replace missing natural source voice in form reed-based element. Behaviour of reed-based element was analysed using computational modelling and using measurements on experimental physical model. The physical model enables changes in setting gap between reed and reed stop and performing acoustical and optical measurements.
|
|
|
Numerical simulation of videokymographic images from the results of the finite element model
Švancara, P. ; Horáček, Jaromír ; Martínek, T. ; Švec, J. G.
The study presents a two-dimensional (2D) finite element (FE) model of the fluid-structure-acoustic interaction during flow induced self-oscillation of the human vocal folds. The FE model combines the FE models of the vocal folds, the trachea and the simplified human vocal tract shaped for phonation of vowel [a:]. The fluid-structure interaction is solved using explicit coupling scheme with separated solvers for structure and fluid domain. The developed FE model comprises large deformations of the vocal-fold tissue, vocal-fold contact, fluid-structure interaction, morphing the fluid mesh according to the vocal-fold motion (Arbitrary Lagrangian-Eulerian approach), solution of unsteady viscous compressible airflow described by the Navier-Stokes equations and airflow separation during the glottis closure. The effect of lamina propria thickness and material properties on simulated videokymographic (VKG) images of vocal-fold vibrations are analyzed. Such variation of the lamina propria properties can be caused by certain vocal-fold pathologies such as Reinke's edema. The developed FE model can be used to study relations among pathological changes in vocal folds tissue, the resulting VKG images and the produced sound spectra.
|
guest ::
